Hidden scale invar iance in molecular van der Waals liquids: A simulation study 



Thomas B. Schr0der, Ulf R. Pedersen, Nicholas Bailey, S0ren Toxvaerd, and Jeppe C. Dyre 
DNRF Centre "Glass and Time, " IMFUFA, Department of Sciences, 
Roskilde University, Postbox 260, DK-4000 Roskilde, Denmark 
(Dated: September 18, 2009) 

Results from molecular dynamics simulations of two viscous molecular model liquids - the Lewis- 
Wahnstrom model of ortho-terphenyl and an asymmetric dumbbell model - are reported. We 
demonstrate that the liquids have a "hidden" approximate scale invariance: Equilibrium potential 
energy fluctuations are accurately described by inverse power law (IPL) potentials, the radial dis- 
tribution functions are accurately reproduced by the IPL's, and the radial distribution functions 
obey the IPL predicted scaling properties to a good approximation. IPL scaling of the dynamics 
also applies - with the scaling exponent predicted by the equilibrium fluctuations. In contrast, the 
equation of state does not obey the IPL scaling. We argue that our results are general for van der 
Waals liquids, but do not apply, e.g., for hydrogen-bonded liquids. 



INTRODUCTION 

A phenomenon is scale invariant if it has no charac- 
teristic length or time. Scale invariance emerged as a 
paradigm in the early 1970's following the tremendous 
successes of the theory of critical phenomena. Liquids 
described by scale-invariant potentials are a theorist's 
dream by having a number of simple properties [l|. If 
potentials vary with distance as r~ n , the configurational 
free energy J 7, as a function of density p = N/V and tem- 
perature T may be written F = NksT f(p 1 /T) where 
7 = n/3. In fact, all properties of inverse power law 
(IPL) liquids - in appropriately defined reduced units - 
are functions of p 1 /T only. Even dynamical properties 
are simple, e.g., the relaxation time r may be written as 

T = t Q g( P yT) (1) 

with the same exponent 7 that characterizes the ther- 
modynamics, to being the reduced time unit (to = 
p^ 1 / 3, y/m/kBT). Unfortunately the predicted equation 
of state does not fit data for real fluids, and IPL liquids 
do not have low-pressure liquid states. Both problems 
derive from the absence of molecular attractions; these 
define the spatial scale of one intermolecular distance at 
low and moderate pressure. Thus real liquids have ap- 
parently little in common with scale-invariant liquids. In 
this paper we demonstrate by example that there is nev- 
ertheless a "hidden" approximate scale invariance in a 
large class of liquids, the van der Waals liquids. Con- 
sequently several properties of IPL liquids apply to van 
der Waals liquids, leading to a number of experimentally 
testable predictions. 

STRONGLY CORRELATING LIQUIDS 

It was recently shown 0, 3, 0] that several model liq- 
uids exhibit strong correlations between the equilibrium 
fluctuations of the potential energy U and the virial 
W (defining the contribution to pressure coming from 



the intermolecular interactions via pV = NUbT + W). 
"Strongly correlating liquids" include [3| the standard 
Lennard-Jones liquid, a liquid with exponential short- 
range repulsion, the Kob- Andersen binary Lennard-Jones 
liquid, a seven-site united-atom model of toluene, the 
three-site Lewis- Wahnstrom model of ortho-terphenyl 
(OTP), as well as a model consisting of asymmetric 
"dumbbell" type molecules. 

For inverse power-law (IPL) potentials oc r~ n there 
is 100% correlation between virial and potential energy 
fluctuations as function of time: AW(t) = jAU(t) where 
7 = n/3 [1]. By reference to the single-component 
Lennard-Jones system it was argued in Ref. [J] that 
strongly correlating liquids are well approximated by in- 
verse power-law potentials as regards their thermal equi- 
librium fluctuations. The present paper presents concrete 
proofs that this is the case - even for molecular model liq- 
uids with van der Waals type interactions (i.e., excluding 
hydrogen, ionic, and covalent bonds). The main conclu- 
sions below are: 

• At a given density and temperature the energy sur- 
face of a strongly correlating molecular liquid may 
be replaced to a good approximation by that of in- 
verse power-law potentials. 

• Strongly correlating liquids "inherit" a number of 
scaling properties from inverse power-law poten- 
tials, but not all such such properties. 

FLUCTUATIONS AND STRUCTURE 

We performed NVT molecular dynamics simulations 
[E 0, 0, S| °f two molecular liquids: i) 512 asymmetric 
dumbbell molecules consisting of pairs of different sized 
Lennard-Jones (LJ) spheres connected by rigid bonds; 
the dumbbells were parameterized to mimic toluene 0]; 
ii) the Lewis- Wahnstrom OTP model 12] (N=324). Both 
models are strongly correlating for the range of state 
points investigated here. For the dumbbell model we 
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FIG. 1: (Color online) (a) AW(t)JN vs. AU(t)/N for 
the asymmetric dumbbell model [9|]. The strong corre- 
lations are quantified by the correlation coefficient, R = 
(AWAU) /^((AW) 2 ) ((AU) 2 } = 0.959. The dashed line 

has the slope 7 = ^ ((AW) 2 ) / ((AU) 2 } = 5.90. (b) 
For configurations taken from the simulation reported in 
Fig- HJa) we evaluated C/ipl.bbW = J2i>j 0iPL,BB(»"ij(t)) 
with ^iPL,BB(r) = Cbb£bb(o-bb/V) 18 . The constant Cbb = 
1.489 (Caa = 1.075, Cab = 1.237) was chosen to make the 
variance of !7ipl,bb(£) equal to the variance of [7bb(£) = 
J2i>j <t>LJ,BB(nj(t)). UipL,BB(t) strongly correlates with 
UBB(t) with a correlation coefficient of 0.977. In compar- 
ison Umff,BB{t) = Ubb (t) — C^ipl,bb(^) has small variance, 
0.045 times the variance of [/bb(£) and is only weakly corre- 
lated with the latter. Redoing the analysis with n = 3 x 5.90 
(the IPL exponent suggested by the fluctuations, see Fig.fjja) 
gives almost identical results. 



find for the correlation coefficient 0.95 < R < 0.97; for 
the OTP model 0.91 < R < 0.92. Fig. [^a) illustrates 
the strong WU-correlation for the asymmetric dumbbell 
model. Since R « 1 it follows that AW{t) « yAU(t) 
in their instantaneous fluctuations 0], consequently 7 is 
referred to as the "slope." 

Fig. [IJb) details the origin of the strong correlations, 
focusing on the strongest interaction in the dumbbell 
model, the 'BB' interaction ('B' indicating the large 
"phenyl" spheres). The fluctuations in C/bb(£) are well 
described by an IPL pair potential <^ipl,bb(?") oc r~ 18 . 
The same conclusion applies for the AA and AB inter- 
actions (data not shown). The reason that <pr>is, q/jM = 
0Lj,a/3(?') — '/'ipl, a/3 (?"*) contributes little to the fluctua- 
tions at constant volume is that this term is approxi- 
mately linear in the region of the first peak of the pair 
correlation function. This implies that, when one pair- 
distance is increased, other pair-distance(s) decrease by 
approximately the same amount, keeping the change in 
U-DiB,af3{t) = J2i>j <k>i8,ap(rij(t)) small Q. Switching 
the ensemble from NVT to NpT (p=1.32GPa) reduces 
the WU-correlation from 0.959 to 0.866, illustrating the 
importance of the constant volume constraint (the im- 
portance increases with decreasing pressure). 

The fluctuations tell us that the relevant part of the 
potential energy surface is well approximated by a cor- 
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FIG. 2: (Color online) (a) Radial distribution functions for 
the asymmetric dumbbell model [j| (full lines), and the corre- 
sponding r~ 18 IPL model (dashed lines, see text), (b) Scaled 
(?ss(r) for three state points; the two state points with the 
same V = p 5 9 /T (black full line and red long dashed line) 
have almost identical (scaled) structure. When increasing p 
at constant V there is a slight shift of the first peak to smaller 
distances, consistent with the existence of the fixed bond. 
Similar scaling is found for the AA and AB interactions (data 
not shown). 



responding 'IPL' system, i.e., the system where the LJ 
pair potentials are replaced by the IPL potentials dis- 
cussed in Fig.[ljb). To test how far this correspondence 
holds, we simulated the IPL system at the same density 
and temperature. Fig. HJa) compares the pair distribu- 
tion functions. The agreement is striking. To the best 
of our knowledge this is the first proof that the structure 
of a molecular liquid is well reproduced by an IPL liq- 
uid. In contrast, the average pressures of the two models 
are quite different (see Fig. [TJb)), because the IPL sys- 
tem does not include contributions to the pressure from 

Consider now changing density and temperature, con- 
serving the parameter T — p 1 /T. For a pure n = 37 
IPL liquid this means that one is effectively studying 
the same system, but with scaled length and time units 
(times scaled by to, lengths by Iq = p -1 / 3 ) d,]: Two state 
points with the same T have 3N-dimensional potential en- 
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FIG. 3: (Color online) (a) Radial distribution functions for 
the LW-OTP model [12] at three different state points. Two 
of the state points (full black lines and red dashed lines) have 
approximately the same F = p 7 ' 9 /T. The third state point 
(full green line) have the same temperature as the first state 
point and same pressure as the same state point, (b) Same 
as in (a), but in reduced units. The two state points with the 
same F (black full line and red long dashed line) have almost 
identical (scaled) structure. Similar to what is seen for the 
asymmetric model (Fig. [2|, there is a slight indication of the 
fixed bond in the first peak of the two state points with the 
same V 



ergy surfaces that are identical in reduced units. Switch- 
ing back to the LJ system amounts to adding the "diff" 
terms of the potentials, which - to a good approximation 
- simply shifts the whole energy surface by an amount 
that only depends on the volume of the state point 
This changes neither structure nor dynamics, and con- 
sequently these properties are predicted to "inherit" the 
scaling properties of the IPL potentials. In contrast, the 
shift of the energy surface does change the absolute val- 
ues of potential energy and pressure (because the shift 
depends on volume), and these properties are not pre- 
dicted to inherit the IPL scaling. 

The IPL scaling of the structure is tested in Fig. Hlb) 
for the asymmetric dumbbell model, and in Fig. [3] for 
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ture for two molecular models 

correlating. For the Lewis-Wahnstrom OTP model P3] the 
zero pressure isobar covers densities in the interval 1.008 
g/cm 3 to f.089 g/cm 3 . For the asymmetric dumbbell model 
[9( the 2.0GPa isobar covers densities in the interval LI 20 
g/cm 3 to 1.217 g/cm 3 . For the three dumbbell isochores the 
pressures at the lowest temperatures are -0.02GPa, 0.67GPa, 
and 1.28GPa respectively. 



the LW-OTP model. For both models the IPL scaling 
works well; there is only a small effect of the fixed bond 
connecting the LJ spheres. Note that average pressure 
and potential energy do not follow the IPL scaling. 



DYNAMICS 

We now turn to the long-time dynamics. To the degree 
that the IPL approximation holds over a range of state 
points, the long-time dynamics is a function of T = p 1 /T 
where 7 is one third of the IPL exponent n, i.e., the 
"density" scaling of Eq. (1) is predicted to apply. This 
type of scaling has in recent years become a well estab- 
lished empirical fact for a large number of viscous liquids 
13l 14 1 (usually ignoring the insignificant state point de- 
pendence of to)- I n particular, density scaling applies 
for van der Waals liquids, but not for liquids with hy- 
drogen bonds or other directional bonds ID, 3, 17, la ]. 
Other forms of scaling have been suggested, but it has 
been argued that the scaling reflects an underlying IPL 
[IQEalliUlllii], and thus that Eq. (1) is the correct 
scaling. 

A potential complication for the IPL explanation of 
density scaling, is that the IPL exponent suggested by 
the fluctuations (n — 37) varies somewhat with state 
point Q [2^|. This is illustrated for the two models in 
Fig.[D Note that, while the suggested exponent is weakly 
state-point dependent for both models, one finds quite 
different exponents for the two models. For the OTP 
model we test density scaling using the slope found by 
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FIG. 5: (Color online) The reduced diffusion coefficient for 
Lewis- Wahnstrom OTP 12]. Open symbols: D* plotted 
without scaling (7 = 0, D* divided by 100). Filled symbols: 
D* scaled according to Eq. (1) with three different scaling ex- 
ponents: 7 = 9.0 (Upper set of curves, D* multiplied by 10), 
7 = 7.9 (Second set of curves), 7 = 7.0 (Third set of curves, 
D* divided by 10). 



averaging over all state points: 7 = 7.9. For the dumbbell 
model we consider two values of the scaling exponent: 
the slope averaged over all state points, 7 = 6.1, and the 
slope averaged over the three data sets with the smallest 
slopes, 7 = 5.9, i.e., the "best" compromise if the p = 
1.006 g/cm 3 isochore is left out of consideration. We 
apply density scaling to the reduced diffusion coefficient, 
D* = (N/V) 1 / 3 (k B T/m)' 1 / 2 D where m is the mass of 
the molecules and D is the diffusion coefficient estimated 
from the long-time mean-square displacement. 

Fig. [5] applies density scaling to the Lewis- Wahnstrom 
OTP model. The scaling works neither with 7 = 9.0 nor 
with 7 = 7.0, whereas it works well with 7 = 7.9 which 
is the exponent we calculated from the fluctuations (Fig. 
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FIG. 6: (Color online) The reduced diffusion coefficient for the 
asymmetric dumbbell model 0. Open symbols: D* plotted 
without scaling (7 = 0, D* divided by 100). Filled symbols: 
D* scaled according to Eq. (1) with four different scaling ex- 
ponents: 7 = 7.0 (Upper set of curves, D* multiplied by 100), 
7 = 6.1 (Second set of curves, D* multiplied by 10), 7 = 5.9 
(Third set of curves), and 7 = 5.0 (Fourth set of curves, D* 
divided by 10). 



0}. In Fig. [6] density scaling is applied to the asymmetric 
dumbbell model. Density scaling works neither with 7 = 
7.0 nor with 7 = 5.0. Comparing the scaling with 7 = 6.1 
to the data without scaling (7 = 0, open symbols), we 
find quite good data collapse - by far most of the density 
dependence is captured by using 7 = 6.1. With 7 = 5.9, 
however, the data collapse is even better for three of the 
data sets, whereas one data set deviates from the master 
curve comprised of these three sets. The latter is the 
isochore p = 1.006 g/cm 3 , i.e., precisely the one that was 
ignored when choosing 7 = 5.9 (Fig.[4|). 

The above results demonstrate that density scaling is a 
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FIG. 7: (Color online) Dynamics of the LW OTP model [Tj] 
in reduced units for the three state point shown in Fig. a) 
Mean -square displacement for the non-central LJ spheres, b) 
Rotational correlation function for the vector perpendicular 
to the plane of the three LJ spheres, using the first Legendre 
polynomial. 
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FIG. 8: (Color online) Results from equilibrium molecular d 
namics simulations of 512 asymmetric dumbbell molecules 
with, respectively, a strong dipole moment (q = ±0.5e, open 
symbols, three isochores), and zero dipole moment (q = Oe, 
filled symbols, three isochores and an isobar), a) Correla- 
tion coefficients, R = (AWAU) / (AW 2 ) (AC/ 2 ), b) The 
'slopes', 7 = y/((AW) 2 ) / ((AU) 2 ). 



consequence of an underlying IPL scale invariance which 
is revealed by studying thermal equilibrium fluctuations 
at a single state point. Obviously the scaling with a sin- 
gle value of 7 can only be perfect for true IPL liquids. 
For "real" liquids a larger region of state points makes 
the scaling less perfect for two reasons (compare Fig. [51 
7 = 6.1 and 5.9): The apparent IPL exponent (n = 3j) 
changes more, and the relative effect of intra-molecular 
interactions that do not obey the IPL scaling (here the 
fixed bonds) is larger. In this connection it is impor- 
tant to note, however, that even small density changes 
are experimentally relevant. Thus changing the density 
of, e.g., phenylphthalein-dimcthylether by just 2% can 
change the relaxation time by more than two ordes of 
magnitude [23| . 

As a consequence of the underlying IPL scale invari- 
ance, state points with same p 1 /T are expected to have 
not only the same scaled structure and reduced diffu- 
sion coefficient as demonstrated above, every aspect of 
the dynamics should be the same in reduced units. Fig- 
ure [7] shows the reduced mean-square displacement and 
rotational correlation coefficient of the LW-OTP model, 
at the three state points of Fig. [3l The agreement be- 
tween the two state points with same p 1 /T (full black line 
and red crosses) is striking. In contrast, the third state 
point (blue dashed line) has the same temperature as the 



first state point (full black line) and the same pressure 
as the second (red crosses), but very different dynamics. 
Clearly, neither pressure nor temperature controls the 
dynamics - to a very good approximation the controlling 
parameter is Y = p 1 /T. 

We argued above that the IPL scaling properties 
demonstrated is a feature only of strongly correlating 
liquids. To illustrate this we modified the asymmetric 
dumbbell model to make it less correlating. This was 
done by adding charges of q = ±0.5e (e being the ele- 
mentary charge) to the LJ spheres, resulting in a dipole 
moment of 7.0D, i.e., almost 4 times stronger than wa- 
ter. The correlation coefficients and slopes are shown 
Fig. [5] (open symbols), compared to the q = version 
of the model (filled symbols). The addition of the large 
dipole moments decreases the correlation, as expected. 
As before we use the averaged slope as the scaling expo- 
nent suggested by the fluctuations, giving 7 = 4.9 for 
q = ±0.5e. The density scaling is applied in Fig. [5] 
Density scaling with 7 = 4.9 works better than either 
7 = 6.0 or 7 = 4.0, but the data collapse is clearly infe- 
rior to what is seen in Fig. [6] and Fig. [5] (note the smaller 
dynamic range in Fig. [9]). This confirms that density 
scaling is indeed a property of strongly correlating liq- 
uids, consistent with the experimental finding that, e.g., 
hydrogen-bonded liquids do not obey density scaling. 
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FIG. 9: (Color online) The reduced diffusion coefficient, 
D* = (N/V) 1/3 (k B T/m)- 1/2 D, for the asymmetric dumbbell 
model with q = ±0.5e scaled according to Eq. (1), with three 
different scaling exponents: 7 = 6.0 (Upper set of curves, D* 
multiplied by 10), 7 = 4.9 (Middle set of curves), and 7 = 4.0 
(Lower set of curves, D* divided by 10). 



CONCLUSIONS 

We have provided direct evidence for an underlying 
IPL scale invariance in model molecular van der Waals 
liquids, leading to three experimentally testable predic- 
tions for these liquids: i) The density scaling exponent 
7 can be identified from the equilibrium fluctuations 
via the fluctuation-dissipation theorem determining the 
frequency-dependent linear thermoviscoelastic response 
functions ^) different state points with same p 1 /T 

- same values of D* and t q - are "isomorphic" , i.e., have 
same scaled structure pB| . iii) Isomorphic state points 
must have the same relaxation spectra. This is consis- 
tent with recent experimental findings for molecular glass 
formers and amorphous polymers (excluding hydrogen 
bonds): Varying pressure and temperature, state points 
with the same relaxation time have the_same shape of 



the dielectric a- relaxation function |26|, |27| . Finally, the 
hidden scale invariance explanation of density scaling ex- 
plains why one observes density scaling of the dynam- 
ics, without the same scaling applying to the pressure 



The hidden scale invariance is a consequence of the 
strong virial/potential energy correlations [1, 0| charac- 
terizing van der Waals liquids, as well as most or all 
metallic liquids. Thus the results from computer simula- 
tions of two model liquids presented here point to a new 
physics of these large classes of liquids, which is simpler 
than previously believed. 

Useful discussions with Kristine Niss are gratefully ac- 
knowledged. This work was supported by the Danish 
National Research Foundation's (DNRF) centre for vis- 
cous liquid dynamics "Glass and Time." 



[2] 
[3] 
[4] 

[5] 

[6] 

[7] 
[8] 

[9] 



[10] 

[11] 
[12] 



[13] 



W. G. Hoover, M. Ross, K. W. Johnson, D. Henderson, 
J. A. Barker, and B. C. Brown, J. Chem. Phys. 52, 4931 
(1970); W. G. Hoover, S. G. Gray, and K. W. Johnson, J. 
Chem. Phys. 55, 1128 (1971); Y. Hiwatari, H. Matsuda, 
T. Ogawa, N. Ogita, and A. Ueda, Prog. Theor. Phys. 
52, 1105 (1974); D. Ben-Amotz and G. J. Stell, J. Chem. 
Phys. 119, 10777 (2003); C. De Michele, F. Sciortino, and 

A. Coniglio, J. Phys.: Condens. Matter 16, L489 (2004). 
U. R. Pedersen, T. Christensen, T. B. Schr0der, and J. 
C. Dyre, Phys. Rev. E 77, 011201 (2008). 

U. R. Pedersen, N. P. Bailey, T. B. Schr0der, and J. C. 
Dyre, Phys. Rev. Lett. 100, 015701 (2008). 
N. P. Bailey, U. R. Pedersen, N. Gnan, T. B. Schr0der, 
and J. C. Dyre, J. Chem. Phys. 129, 184507 (2008); J. 
Chem. Phys. 129, 184508 (2008). 

NVT simulations were carried out using Gromacs soft- 
ware using the Nose-Hoover thermostat 0. Molecules 
were kept rigid using the LINCS 8] algorithm. 
H. J. C. Berendsen, D. van der Spoel, and R. van Drunen, 
Comp. Phys. Comm. 91, 43 (1995); E. Lindahl, B. Hess, 
and D. van der Spoel, J. Mol. Mod. 7, 306 (2001). 
S. A. Nose, Mol. Phys. 52, 255 (1984); W. G. Hoover, 
Phys. Rev. A 31, 1695 (1985). 

B. Hess, H. Bekker, H. J. C. Berendsen, and J. G. E. M. 
Fraaije, J. Comp. Chem. 18, 1463 (1997). 

The large LJ sphere, mimicking the phenyl group, is sim- 
ilar to the one in the Lewis- Wahnstrom OTP model (TTJ] 
with the parameters m p — 77.106 u, a v — 0.4963 nm 
and e p = 5.726 kj/mol. The small sphere, mimicking 
the methyl group, was taken from UA-OPLS [l"H with 
m m = 15.035 u, a m = 0.3910 nm and e m = 0.66944 
kj/mol. Bond length: d = 0.29 nm. The interaction 
between unlike particles is determined by the Lorentz- 
Berthelot mixing rules. 

G. Wahnstrom and L. J. Lewis, Physica A 201, 150 
(1993); L. J. Lewis and G. Wahnstrom, Phys. Rev. E 
50, 3865 (1994). 

W. L. Jorgensen, J. D. Madura, and C. J. Swenson, J. 
Am. Chem. Soc. 106, 6638 (1984). 

N — 324 molecules consisting of three LJ particles (with 
a = 0.483 nm, e = 600 Kk B ~ 4.989 kj/mol and m = 
76.768 u) placed in the corners of a rigid isosceles triangle 
with two sides of length a — 0.483 nm and one angle of 
75° [lOj . LJ potentials were cut at r c = 2.5<r using a shift 
function [6j] with ri = 2.3a. 

A. Tolle, Rep. Prog. Phys. 64, 1473 (2001); C. Dreyfus, 
A. Aouadi, J. Gapinski, M. Matos-Lopes, W. Steffen, 



7 



A. Patkowski, and R. M. Pick, Phys. Rev. E 68, 011204 

(2003) ; C. Alba-Simionesco, A. Cailliaux, A. Alegria, and 
G. Tarjus, Europhys. Lett. 68, 58 (2004); C. M. Roland, 
S. Hensel-Bielowka, M. Paluch, and R. Casalini, Rep. 
Prog. Phys. 68, 1405 (2005). 

[14] R. Casalini and C. M. Roland, Phys. Rev. E 69, 062501 

(2004) . 

[15] A. Grzybowski, K. Grzybowska, J. Ziolo, and M. Paluch, 

Phys. Rev E 74, 041503 (2006). 
[16] C. M. Roland, S. Bair, and R. Casalini, J. Chem. Phys. 

125, 124508 (2006). 
[17] A. Le Grand, C. Dreyfus, C. Bousquet, and R. M. Pick, 

Phys. Rev E 75, 061203 (2007). 
[18] C. M. Roland, R. Casalini, R. Bergman, and J. Mattsson, 

Phys. Rev. B 77, 012201 (2008). 
[19] C. M. Roland and R. Casalini, J. Phys.: Condens. Matter 

19, 205118 (2007). 
[20] N. P. Bailey, T. Christensen, B. Jakobsen, K. Niss, N. B. 

Olsen, U. R. Pedersen, T. B. Schr0der, and J. C. Dyre, 

J. Phys.: Condens. Matter 20 244113 (2008). 
[21] D. Coslovich and C. M. Roland, J. Phys. Chem. B 112, 

1329 (2008); J. Chem. Phys. 130, 014508 (2009). 



[22] T. B. Schr0der, U. R. Pedersen, and J. C. Dyre, 

larXiv:0803.2"T99 (2008). 
[23] M. Paluch, R. Casalini, and C. M. Roland, Phys. Rev. B 

66, 092202 (2002). 
[24] N. L. Ellegaard, T. Christensen, P. V, Christiansen, N. 

B. Olsen, U. R. Pedersen, T. B. Schr0der, and J. C. Dyre, 

J. Chem. Phys. 126, 074502 (2007). 
[25] N. Gnan, T. B. Schr0der, U. R. Pedersen, N. P. Bailey, 

and J. C. Dvre. larXiv:0905. 34971 (2009). 
[26] K. L. Ngai, R. Casalini, S. Capaccioli, M. Paluch, and C. 

M. Roland, J. Phys. Chem. B 109, 17356 (2005). 
[27] C. M. Roland, Soft Matter 4, 2316 (2008). 
[28] A. Grzybowski, M. Paluch, and K. Grzybowska, J. Phys. 

Chem. B 113, 7419 (2009). 
[29] Another issue is: what is the optimal exponent to use 

when correlation is not perfect? Any exponent ■yR k (k 

being a constant) has the required property of reducing 

to the pure IPL exponent in the limit of R going to one. In 

the present work we chose to test the exponent 7 defined 

in Fig. [j] 



